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ABSTRACT 

A thermally induced spin rate ripp*e hypothesis is proposed to explain the spin rate anomaly 
observed on ISEE-B (launched October 22, 1977). This involves the two radial 14.5 meter beryllium 
copper tape ribbons going in and out of the spacecraft hub shadow. A thermal lag time constant is 
applied to the thermally induced ribbon displacements which perturb the spin rate. It is inferred 
that the averaged thermally induced ribbon displacements are coupled to the ribbon angular 
motion. Qualitative analytic results show a possible exponential build up of the inplane motion of 
the ribbon which in turn causes the spin rate ripple, ultimately limited by damping in the ribbon 
and spacecraft. Qualitative increase in the oscillation period is also indicated and that the thermal 
lag is fundamental for the period increase. Numerical parameter values required to agree with 
in-orbit initial exponential build up (after a torquing maneuver) are found to be reasonable; those 
required for the ripple period are somewhat extreme. 
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THERMALLY INDUCED SPIN RATE RIPPLE 
ON SPACECRAFT WITH LONG RADIAL APPENDAGES 

I. INTRODUCTION 

On October 22, 1977 the ISEE-B spacecraft, which is part of an international program between 
NASA and ESA, was launched into a highly eccentric earth orbit. After all spacecraft appendages 
were deployed and time had elapsed to damp out deployment transients, a disconcerting phenome- 
non was noted in the spin rate. There appeared to be sustained spin period fluctuations. Attempts 
were made to eliminate the spin rate ripple with on-board thruster and ground timed commands. 
Although the spin period fluctuations were significantly reduced by the torquing maneuver, it 
returned after a period of time to full pre-torquing magnitude. Another puzzling aspect of the spin 
rate ripple was the period. Instead of the expected antisymmetric mode period of 12.8 seconds 
associated with two radial 14.5 meter tape ribbon booms, a period of 13.2 seconds was observed. 
The symmetric mode vibration period, which was not observable with the on-board sun sensor, had 
a period of 13.5 seconds. Longman and Massart in Reference l give a resume of events that took 
place. After examining possible perturbing effects, they felt that solar pressure acting on the two 
radial 1 4.5 meter tape booms having small tip masses might be the cause of the spin rate ripple. 
Subsequent theoretical analysis and computer simulations did not sustain the solar pressure effect 
hypothesis. It is the intent of this paper to show that thermal contraction and expansion of the 
14.5 meter ribbon booms going in and out of the spacecraft hub shadow could be a reasonable 
explanation of the observed spin rate ripple phenomenon. The sun’s rays were essentially perpen- 
dicular to the spin axis of the spacecraft, at times perpendicular to the wide surface of the ribbon 
during a rotation period and hub shadowing of the ribbons does take place. The equations of 
motion and related thermal equations for this situation will now be developed. 
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2. EQUATIONS OF MOTION 

As stated in Reference 1 , the observed sun sensor data from ISEE-B indicated that there was no 
steady state nutation motion, that the dynamic motion of the spacecraft and flexible booms was 
essentially planar. This greatly simplifies the mathematical formulation and analysis of the problem. 
Figure 1 shows a sketch of the spacecraft with appendages. 



The main focus of attention is the rigid hub and the two 14.5 meter ribbon booms. The ribbon 
booms are made of beryllium copper with a width of 0.5 centimeter (0.1 97 inch) and a thickness 
of 0.04 centimeter <0.0157 inch). Figure 2 defines the planar coordinate systems used in this 
analysis. 



ATTACHED TO HUB 
NON ROTATING 


INERTIALLY 

FIXED 


Since a force-free system is being considered, the origin of the reference coordinate system is 
located at the center of mass of the hub-booms system. The usual assumption of attitude motion 
uncoupled from the spacecraft orbit is made and so the center of mass is considered fixed (or 
moving with constant velocity). The spacecraft hub is allowed two linear degrees of freedom in a 
plane plus rotation. The ribbons and tip masses are modelled as physical pendula hinged at their 
base. This latter consideration models the lowest vibration modes of the ribbon with reasonable 
accuracy. 

The Lagrangian method will be used to obtain the equations of motion. This requires determining, 
in this case, only the total kinetic energy of the system. Let R be a vector from the center of mass 
of the system to any point on the spacecraft and R the total time derivative of the position vector, 

m 

then the total kinetic energy of the system and also the Lagrangian is 


3 


V 
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T *“* / R • R dm ( 1 ) 

2 Jsic 

where dm is an increment of mass and the integration is over the entire spacecraft. To simplify the 
analysis, in the derivation of the kinetic energy of the spacecraft, the contribution of fi to the total 
kinetic energy will be neglected, but 9 in the equations of motion will be retained. It will be seen 
that this is adequate to reveal the phenomenon under consideration. Let R 0 be a vector from the 
center of mass to the center of the hub and 7a vector from the hub center to a point on the space- 
craft. Then the following can be written 

R = R 0+ ? 

Rq a W + yj, R 0 3 W + fj 

7 s * xi + yj,7= xi + yj 

hence 

R = (77 + x)i + (r + y)j (2) 

where i and j are unit vectors as shown in Figure 2. In taking the time derivatives, there is no u 
cross terms because the coordinate systems are not rotating. Substituting Equation 2 into Equa- 
tion 1 and canying out the vector dot product results in 

T / (fr+f 2 )dm +-;•/ (x : + y 2 )dm + / (t?x + fy)dm (3) 

--/s/c -Jsic Jsic 

It will be noticed in Equation 3 that the first integral relates to translation, the second integral 
relates to rotation and the third integral relates to translation-rotation coupling. 

Since 77 and £ are independent of the spacecraft mass distribution, the first integral can be integrated 
to give 

Xc (T > 2 + ^ )dm = m t ^ 2 + (4) 

where Mj is the total mass of the spacecraft. 


4 
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In general, the dynamic system that we are considering has the degrees of freedom: 77 , a, x and y 
or as will be seen subsequently, 77 , a, and . The third integral in Equation 3 can be evalu- 
ated in the following way: the equation of motion for the 77 direction using the Lagrangian method 

is 


d / 3T \ dT 




dt \ 3?7 / 877 7 


Now from Equation 3,—* * 0, hence Equation 5 can l* integrated to 

3*7 


3T 

— = constant 
dr) 


3T 

Evaluating — from Equation 3 and 4 gives 
377 


3T 

377 



dm = constant 


(5) 


( 6 ) 


(7) 


At t * 0, everything is assumed quiescent: 77 3 x = 0, hence the constant in Equation 7 is found to 
be zero and 




~ f SIC 


x dm 


( 8 ) 


In a similar manner for the f linear direction results in 


fM T*‘Xc i ' dm 


(9) 


Using Equations 8 and 9 and noting that 77 and £ do not depend on the integration over the S/C, the 
third integral in Equation 3 evaluates to 


f (t?x + fy)dm = - M t (t7 2 + f 2 ) 
J sic 1 


( 10 ) 


The Lagrangian thus far can be written as 
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Continuing further to eventually express rj and f in terms of angular displacements of the ribbon so 
that the Lagrangian can be completely expressed in terms of these angular displacements. Equations 
3 and 9 can be integrated once more to give 


and 


” M T "L 


x dm 


(12a) 



(12b) 


Note that Equations 1 2a. b are simply the first mass moment equations. Using the transformation 
Xj * a cos a + $ cos (0j + a),x n ■ a cos (a + n) + £ cos (6-, + a + ir) 
yj * a sin a ♦ £ sin (5 1 + a), y 2 * a sin (o + ir) + £ sin (0 2 + ck + tt) (13) 


where a is the radius of the hub, £ is the distance along the wire to the increment of mass and 
canying out the integration indicated by Equations 1 2a and 1 2b (the hub does not contribute): 



wire & 
tip mass 



px 2 d$ 


wire & 
tip mass 


(13a) 


fM T py^S-J py : d£ 

wire & wire & 

tip mass tip mass 

where p is the nominal mass per unit length of the wire, results in 


(13b) 


*- (p£ 2 / 2 + m£)[cos(0j +a)- cos (0-, + a)] (14a) 

fM-j- * - (p£ 2 /2 ♦ m£)(sin (0j + a) - sin (0 2 + Q )1 (14b) 

where an average wire length, £, has been used to simplify the algebra for this part of the kinetic 
energy and m is the tip mass. 



ORIGINAL F4GE 19 
OF POOR QUALITY 

Differentiating Equations 14a and 14b and dividing by Mj give} the spacecraft translation velocities 
in tenns of the angular motion of the hub and ribbons 


jj • ^r, [($ + a) sin (8 , + a) - ($■, ♦ a) sin (8? + a)J 

1 * x 

(o2 2 /2 + m2) 

t ■ - ■ ((5 j + a) cos (8 j + or) - ( 8 2 + a) cos (0 2 ♦<*)) 

•' 1 ^ 

Squaring Equations 15a and 15b, adding and multiplying by M T gives 


(15a) 


(15b) 


M t (t»- + f z )«M c [(0, + d) 2 + ( 8 2 + d)- - 2(0 , + d)(0 2 +«)] cos (0, - 9 2 ) (16) 


where M c is defined as 


?j + 

and 2 is an average length, — " ■ ~ 


M. 


(pg-/2 * m2)* 
M t 


(17) 


In a somewhat similar manner, the integral in Equation 1 1 related to the rotational kinetic energy, 
T R , can be evaluated 

T r »-i-Idr +~[(p2j + mj)a-6r + 2(pS.]/~ + mj)acir(6j + d) cos 8 l 
♦ (mj2j /3)(8 { + d) 2 ] 

l(p« 2 + m 2 )a 2 or + 2 (pil /2 + m2-,) ad(tf 2 + «) cos 8-, 

+ (m 2 2^+,o2;/3)(fl 2 +i) 2 ] (lS 

In Equation 18, 1 is the hub moment of inertia and an attempt has been made to include assymme- 
tries caused by different ribbon lengths and tip masses. Initially it was thought that the small dif- 
ference in tip mass weight might be a contributing factor to the anomaly. This is briefly discussed 
in Section 5. 


Defining the following quantities 


7 
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a 2 (pfij + mj) 

M 2 j * a 2 (p £ 2 * m 2 ^ 


afpfijC + m^j) 

M 22 ■ a (pi\/2 + m 2 £ 2 ) 


pfif/3 + mjfif 

M 23 a pill 3 + on £2 

(19) 


and using Equations 1 1, 16, 18 and 19, the Lagrangian (total kinetic energy) can be written as 

T *--I 6 r +-^[Mjjor + ZM^cWdi + a) cos dj + M 13 (d, + a) 2 ] 

«• * 

+-“[M2iCr + 2M22<*(02 * cos & 2 + ^23^1 + 

«• 

- --M c C(^ i + a)' + (d-, + a ) 2 - 2 (dj + a )(d, + a) cos ( 6 j - d-Ol 
Tlie remaining angular equations of motion can now be written. For 9 , , 


CO) 


d / 3T V 3T 

,-!♦ I za ) $0 * 


dt V30J 


* -cd , 
1 1 


Cl a) 


or in integral form 


ii.iil ♦ n 

39 1 39 ,J J \ 

t-0 


st ;-'- 9 * 1 dt 


(2 lb) 


where 


3T 

3d, 


Mj 2 « cos dj + M 13 (dj + a) - M c [dj + a- (6-, + a) cos (d, - d-,)l 


(22) 


3T 

3d, 


'-Mj ->(dj + a) sin dj - M c (dj + a)(d, + a) sin (dj - d-,) 


(23) 


and-cdj is a damping term due to combined ribbon and spacecraft damping. For d-, , 


or in integral form 


d / 3T \ 3T . 

T[le7! -SdT'V'^ 


3T 3T1 r x /b t . \ J 

dt 


(24a) 


(24b) 


t-0 
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where 



■ co» 9 2 + + <*) " M c (02 + *“ (^j + a) cos <6 j - 02 ^ 


(25) 


3T 


— — *. + a> sindi -t* i + d)(0'» + a) sin ( 8 , - 9 -,) 
30 , .ci . i - 


and -c0-> is a boom motion damping term. 


(26) 


For a. 


where 


and 


d / 3T \ 3T 
dt \ 3d / 3a 


Q a -0 


-jt* la + Mj j<i + Yj] ( 8 1 + 2a) cos 9 j + M| ^(^] + a> 

+ M-, ja+ MWtN + 2a) cos 9 ^ + ^ 23(^2 + ° t) 

- M c l(0j + d) + (82 + a) - (0j + 9 z + 2d) cos (9 l - 9 2 J 


( 27 ) 


(23) 


3T 

3a 


j. 


Now Equation 27 can be integrated directly because of Equation 29 ; this results in 


(29) 


Id + M j j a ♦ M j j (9 j + 2d) cos 8 j + M 2 3 (0 { + a) 

+ M-» jd + + 2d) COS 9 2 + M2 3 (^2 + ®) 

- M c [(0j + a) + (0 ; +a)- (5j +9 2 + 2d) cos(0j - 0 ; 1 * constant ■ H 0 # (30) 
which is the conservation of angular momentum equation, since no external torques are considered 


acting. 
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3. linearized equations of motion 

The partially linearized equations of motion (derived from Equation 21b and 24b) taking into 
account time variations of ribbon length can be written as 

Mj 3(0*1 +a) + Mi2« + (Mi3 •+• M j 2 “ M c (^i " $2)“ ^ 

+ c0j + Mi 2 ® £ *'®i + M c cr(0 . - d - i ) = 0 (31) 

and 

M 2 3 (0*2 + 5) + M 22 « + (M 23 + M 22 )a + M c (0, - S 2 ) + M c ( 0 j - B 2 ) 

+ c0 2 + M2 2 a 2 02' M c“ 2 ^l '®2 ) = 0 (32) 

From the definition of the M’s, we note that 

Mi3+Mi2=’i(lMi2 + Mn)gi (33a) 

+ M-»-» =- (2M^ + Mo 1 )« 2 . (33b) 

It is convenient to express the linearized equations of motion in terms of normal modes of vibration: 
antisymmetric mode. 0j = 9 j + 0 2 > ^ symmetric mode, t - B-, . Adding Equations 3 1 and 
32 and considering 


M-) 3 ^ M 1 3 , M 2 2 588 M j 2 j 


both constant and from the linearized conservation of angular momentum (Equation 30) 

l s x s 

we can express the 0j mode differential equation as 

... 2 1 • 

I N 0! + c0j + M, 2 0J 0 01 = ( 2 Mj 2 + Mj jKSj + 

3 


(34) 


(35) 


10 


where 
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*c 


and cj 0 is the nominal spin rate. For the <t>-, mode differential equation we can write 


(M| 3 - + c^2 + ^1 Z + 2 M C ) w o“*2 = (^Mj2 + Mj j )(2j - 

*** ci 


(36) 


Both equations axe considered to have constant coefficients. From the right hand side of Equations 
35 and 36 we see the potentially driving nature of time varying ribbon lengths that can come from 
contraction and expansion of the ribbon booms going into and out of the spacecraft hub shadow. 


4, THERMAL EFFECTS 

Additional equations are required for the contraction and expansion of the ribbon going in and out 
of the shadow of the spacecraft hub. Initially the ribbon will be divided up into small typical ele- 
ments and shadowing effects brought in by averaging over a spin period: subsequently the elements 
are sumed for the total thermal displacement. The thermal effect will be modeled as a first order 
lag differential equation of the form 

i s 

5,+-^=— Fj (0) (37) 

J r J r i 

where 5j is the thermal change in length of a small ribbon element, r is the characteristic time con- 
stant of the thermal effect, 5 C is the steady state value and F ; is a shadowing factor: 1 in the sun- 

s j J 

iight and 0 in complete shadow; 0 is the angular location of the small element about the spin axis 
(see Figure 3). For simplicity, conduction along the ribbon is neglected. Equation 37 is related 
to a linearization of a heat balance between heat storage loss and radiation to space plus the rela- 
tionship between thermal displacement and temperature difference. 



We are mainly interested in long term effects. Equation 37 will be averaged over a hub rotation to 
take out the time consuming spin rate dependence. This procedure essentially retains the sym- 
metric and antisymmetric frequencies (which are lower than the spin rate by a factor greater than 
four) and effectively results in using only the “constant” term in a Fourier series expansion of 
Fj 03). Averaging Equation 37 results in 


i + L S « JL ^F 03)da 

J r J r2?r J 1 


(31 


Referring to Figure 3 the integral on the right hand side of Equation 38 can be written as 


i / i j r* r* r 2 * i ( 

-r- 0Fj(/3)da = — - \ / Ida + / Oda + I 1 da = — - < a. -a + lit 

2?r jT J 2w (Jo J 2ir\ s i s e 


(3 
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where a s . is the hub rotation angle going into the shadow of the j tb element and is the exiting 
angle from the shadow. Now the onset of shadowing for each element is given by the following 
condition, linearized for small ribbon angles 

(a + £ s ) cos 0 S . = a 

or 

a 

3 arcos (40) 

s i a+| s 

where £ s identifies the element going into the shadow. The linearized condition for exiting the 
shadow is given by 


(a + | ) cos = -a 
s 5 e 


or 


( 3 . = arcos 
s e 



For small ribbon angles, the following geometric relations can be written from Figure 3 


„ «.«i 


or 





(41) 


(42) 


and 


or 




s e 


l 3 a + T e * a + “ 

s e s e e s i a + 





3 + *s 


( 43 ) 
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Using Equations 40, 41, 42 and 43 and forming the difference aq - a e , we arrive at 

u 


a \ /a 

a - a. = arc os ( — - arcos - 

s i s e \ a + L J \ a • 


*TJ * a+T, (# * " $>> 


(44) 


Since 


arcos (x) + arcos (~\)~n 


(45) 


the above equation can be written 


<*, - a. 


arcos 




— + (d e -9:)-n 

+ L a + L 


(46) 


•s ■' ~ S s 

Now strictly speaking, 9 i and d s vary somewhat as one goes along the ribbon length entering and 
exiting the shadow. Neglecting this effect and defining as the net angular motion gen- 

erated going through the hub shadow, Equation 38 (using Equations 39 and 46) for ribbon 1 
becomes 


5 S; 


«lj + 7«U*2^1' + 2a,COS 


- J£l| 


a + / a+ 


(47) 


For ribbon 2 we have 


5-i; + — 5-i; = 


s j 


-J 


2J 


lltT 


rr + 2 arcos 


a + £< 




a + £< 


(48) 


Basically the right hand side of Equations 47 and 48 (braces divided by 2ir) is a measure of the frac- 
tion of a rotation period that is spent in the sunlight by an element. The minus sign is chosen for 
the net angular motion (and also guided by the subsequent results) because as the ribbon enters the 
shadow, Coriolis forces (caused by contraction) urge the ribbon forward into the shadow; as the rib- 
bon comes out of the shadow, Coriolis forces (now caused by expansion) tend to keep the ribbon in 
the shadow; this effectively reduces the time in the sunlight and the negative sign reflects this. Thus 
we see that the ribbon angular motion can affect the thermally induced displacements by the time 
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spent in the sunlight/ shadow. This is a possible mechanism for solar thermal energy to be delivered 
to the dynamic system. 


Summing over all of the small elements to obtain the total ribbon thermal displacements gives 




V = ^ 2 


(49a) 

(49b) 


Note that 2j = 2 Q + 5j and i 2 * 2 0 + 5 2 where is the nominal length of the ribbon. Equations 
47 and 48 become 

M 


| ff + 2arcos (lT^) T7 

1 1 _ ^ 

5^ +— , 5-» a 5 : 

- r * 2ttt L~d SJ 


7T + 2 arcos 


/ a+ £. 


(50a) 

(50b) 


Now the right hand side of Equations 50a and b can be converted to integrals by noting that 


5=5 

Si U 5 i £ 


*J 


(51) 




Each infinitesimal element contributes a fraction,——, to the total steady state change of length, 5 . 


We thus have 


and 


1 5 s ( 2 r / a \ . 1 f'Wl \ 

5, +-*-6, 3 r * + T / arcos ( - ) d$-- / df 

1 t 1 2tt ( tJ 0 \a + |/ i J 0 a + £ s 

• 1 . 5 s | 2 /• / a \ I p .J 

5-, +— 5, = - — W+’r/ arcos/ H df--/ T d W 

’ r 2 2 tt ( 2./ \a + $/ V 0 a + * s j 


(52a) 


(52b) 


where 


Z 5 Si =5 S 


(53) 
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the total steady state change of length of the ribbon, and integration is indicated for the second and 
third terms. Carrying out the integrations and defining quantities results in the compact form 


5, +—5, *-—(e- k0, ) 
1 t 1 t 1 


§2 +-“6 2 3 ( e ~ k0 2 ) 


(54a) 

(54b) 


where 


e * 


£{• 


rr + 2( 1 + a/2) arcos 


a 2a i (a + 2 + ^/i‘- + 2aJ i ) ^ 


a + 2 2 


) 


(55) 


and 


k= T“ I--~ln(l+2/a)i 
2n l Z \ 

The length 2 considered here is the nominal length of the ribbon. Using operator notation 
(d/dt a D), Equations 54a, b can be written as 


(56) 


/ 1 \ 

5 S 


( d+ t) 

5j (e- k6j) 

(57a) 


5 S 


( d+ 7) 

5* * — (e - k0->) 

i T 

(57b) 


and the component of the averaged steady state solution depending on the ribbon angles can be 
expressed symbolically as 


6 1 


5. ■ - 


5 s k *l 
(rD+ I) 

5 s k0 2 
(rD+ 1) 


(58a) 

(58b) 


We can thus write the sum and difference of the thermal effect of the ribbons in terms of the , 
<t>2 modes 


5j + 5 2 = + 2 2 


5 s k< *l 
rD + 1 


(59a) 
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and 


^ 1 ~ ^2 * 1 “ ®2 


5jk02 
rD + 1 


(59b) 


5. ANALYSIS AND INTERPRETATION OF RESULTS 

We are now in a position to return to the linearized equations of motion and evaluate the effect of 
the thermally induced ribbon displacements. Substituting Equation 59a into Equation 35 and 
Equation 59b into Equation 36 results in the governing dynamic-thermal differential equations for 
the 0 j and 0-, modes 


! n 0i 


and 


(Mjj - 2M.)0 2 + 


5 s kw 0 (2M i: +M n ) 
a (rD + 1) 


5 s k(*> 0 (2Mj 2 + Mj j) 


0j + = 0 


(60a) 


0-, + (Mj2 + 2M c )cj o "02 = 0 (60b) 


(rD + 1) 

Focusing on the antisymmetric mode, 0 { , the nature of the motion can be ascertained by assuming 
r small and expanding the operator 


1 


1 + rD 


• = 1 - rD . . 


(61) 


in Equation 60 a. After rearranging one can write 
— ku; 0 (2Mj 2 +M U )- c 

cl 


0!-' 


M I2 w 0 ^1 


r N + 


rw 0 5 s k 


’l 


(2M 12 + M u ) 


h + 


rw 0 5 s k 


(62) 


(2M 12 + M h ) 


Let 


■ko; 0 (2M 12 + M n j- c 




In + 


rw o 5 s k 


(2Mj 2 + Mj j) 


(63) 
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and 


& 


M 12 w o‘ 


In- 


rw 0 5 s k 


(2M 12 +M n ) 


(64) 


Equation 62 can be put into the standard 2nd order differential equation 


0j - 3 0 


(65) 


which has the solution 


0i »Ae^ flt sin i/l - M 2 ftt + Be^ nt cos 1 - 


( 66 ) 


where A and B are constants. For combined ribbon and spacecraft damping less than the thermal 
driving term 

c < 5 s kw 0 (2Mj 2 +M, j)/a 

M is positive, we note that the amplitude of the antisymmetric mode increases exponentially, also 
note from Equation 64 that the frequency is reduced by the r term. A similar result can be 
obtained for the symmetric 0 2 mode. 


Figure 4 shows some spin ripple amplitude flight data (taken from reference 1) after a torquing 
maneuver on April 23, 1980. Immediately after the maneuver, a divergent time constant of approx- 
imately 3.6 days was calculated. To show some quantitative results from this analysis, Equation 
60a was cleared of the 1 + rD operator in the denominator of the first derivative tenn and only 
linear terms were retained. This results in a third order differential equation for 0j , 

rI N 0j + (I N + ct)0j + (c + rMj 2 cj q 2 - G)0j + Mj 2 cu o 2 0j = 0 (67) 


where 


SPIN RIPPLE 
AMPLITUDE 
(COUNTS) 


* 

* 
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I wun; 4 . liyolulio,, ofSpin Ripple Amplilude follow, eg Ripple Rcrluclion MnncimM Reference I) 
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For the following spacecraft parameters 

I » 62.033 slug ft“ p * 2.8 X 10” 5 slug/ft 

I s * 69.79 slug ft 2 5 S » 0.15 ft 

a » 2.082 ft u> 0 ■ 2.0724 rad/sec 

e 0 » 47.54 ft lc* 0.159 

c » 0.0 

the roots of the ^ mode are displayed m Table I as a function of the thermal time constant, r. 


Table I 

Roots of Thermally Induced Spin Rate Ripple Equation - Antisymmetric Mode 


Tau 

Root 1 

Root 

■) 

Root 3 

Real 

Imag. 

Real 

Imag. 

Real 

Imag 

.0 

2.5477E-07 

HRI 

2.5477E-07 

-.19120 

-6.6672E-03 

.0 

5.00000 

1.9565E-04 

1 

1.9565E-04 


-.20039 

.0 

10.0000 

5.5097E-05 


5.5097E-05 

-.49095 

-.10011 

.0 

15.0000 

2.4842E-05 

■ 

2.4842E-05 

-.49104 

-6.6716E-02 

.0 

20.0000 

1.4097E-Q5 

.49108 

1.4097E-05 

-.49108 

-5.0028E-02 

.0 

25.0000 

9.0796E-06 

.49111 

9.0796E-06 

-.49111 

-4.0018E-02 

.0 

30.0000 

6.3078E-06 

.49113 

6.3078E-06 

-.49113 

-3.3346E-02 

.0 

35.0000 

4.6525E-06 

.49114 

4.6525E-06 

-.49114 

-2.8581E-02 

.0 

40.0000 

3.5606E-06 

.49115 

3.5606E-06 

-.49115 

-2.5007E-02 

.0 

45.0000 

2.8137E-06 

.49116 

2.8137E-06 

-.49116 

-2.2228E-02 

.0 

500000 

2.2306E-06 

.49117 

2.2806E-06 

-.49117 

-2.0005E-02 

.0 

55.0000 

1.8836E-G6 

.49117 

1.8836E-06 

-.49117 

- 1.8186E-02 

.0 

60.0000 

1.5833E-06 

.49118 

1.5833E-06 

-.49118 

-1.6670E-02 

.0 

65.0000 

1.3459E-06 

.49118 

1.3459E-06 

-.49118 

- 1.5387E-02 

.0 

70.0000 

1.1642E-06 

.49118 

1 . 1 642E-06 

-.49118 

- 1 4288E-02 

.0 

75.0000 

1.0129E-06 

.49118 

1.0129E-06 

-.49118 

- 1.3335E-02 

.0 

80.0000 

8.8824E-07 

.49119 

8.8824E-07 

-’49119 

- 1.2502E-02 

.0 

85.0000 

7.8724E-07 

.49119 

7.8724E-07 

-.49119 

- 1.1766E-02 

.0 

90.0000 

7.025 7E -07 

.49119 

7.0257E-07 

-.49119 

-1.11 13E-02 

.0 

95.0000 

6.2844E-07 

.49119 

6.2844E-07 

-.491 19 

- 1.0528E-02 

.0 

100.000 

5.6786E-07 

.49119 

5.6786E-07 

-.49119 

-1.0001E-02 

.0 

105.000 

5.1 158E-07 

.49120 

5.1 158E-07 

-.49120 

-9.5248E-03 

.0 

110.000 

4.7019E-07 

.49120 

4.7019E-07 

-.49120 

-9.0919E-03 

.0 

115.000 

4.2973E-07 

.49120 

4.2973E-07 

-.49120 

-8.6965E-03 

.0 

120.000 

3.9520E-07 

.49120 

3.9520E-07 

-.49120 

-8.3341E-03 

.0 

125.000 

3.6588E-07 

.49120 

3.6588E-07 

-.49120 

-8.0007E-03 

.0 

130.000 

3.3883E-07 

.49120 

3.3883E-07 

-.49120 

-7.6930E-03 

.0 

135.000 

3.1263E-07 

.49120 

3.1263E-07 

-.49120 

-7.4080E-03 

.0 

140.000 

2.905 IE-07 

.49120 

2.905 IE-07 

-.49120 

-7. 1434c -03 

.0 

145.000 

2.7096E-07 

.49120 

2.7096E-07 

-.49120 

-6.8971E-03 

.0 

150.000 

2.5477E-07 

.49120 

2.5477E-07 

-.49120 

-6.6672E-03 

.0 


ooooooooooooooooooooooooo 
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It win be noted thn for a r of 40 seconds (which does not appear to be unreasonable) a divergent 
time constant of 3.3 days is possible. In regard to the frequency of oscillation, the frequency is 
lower (period longer) than without thermal excitation, but it is not is low as what was observed in 
orbit. 

Referring back to Figure 4, it will be noted that the rate of spin ripple amplitude growth decreases 
with time and eventually stops growing; an essentially “steady state” condition occurs. Now it is 
inferred that the damping parameter, c, is slowly varying, increasing as the amplitude of ribbon 
motion increases, eventually limiting the motion. A quasi steady state result can be obtained from 
this analysis that clearly shows the increase in the antisymmetric mode period. The condition is 
imposed that the roots of liquation 67 have two purely imaginary roots and one real. This implies 
a “steady state” sinusoidal motion after the transient damps out. The cubic equation for the roots 
of Equation 67 should have the form 

(s + r t )(s~ + ft“) * r* + rj s- + £2*s + r, * 0 (69) 


The conditions are imposed that 

1 c 


r, *-+r- • «— + 


, c w o" M 12 G 


r I 


N 


rI N 


*N 


and 


. m 12 w o‘ 

r ,r.2.— - 


Eliminating c from the equations and solving for the frequency results in a quadratic in £2“ : 


r 2 ^. 4 + (1 + Gt/I n - r^cjj 2 )^ 2 - ■ 0 


(70) 


where 


Wi 


2 M 12 w c 
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the natural frequency of the antisymmetric mode without the thermal effect. Letting 

R*£2 2 /wf (7i) 

Equation 70 can be put into the form 

wj^nR 2 + ^ 1 +“ — rw]‘ j R - 1*0 (72) 

A simple approximate expression can He obtained for R which reveals the nature of the thermal 
effect on the natural frequency of the antisymmetric mode by considering 


Rvl+e 


(73) 


where e is much less than one. Substituting Equation 73 into 72 and neglecting e- and solving for 
e results in 

Gr/I N 


e *- 


1 + Gr/l y + uij'r 2 


(74) 


Thus 


£2* w cj,* 1 - 


Gr/U 


^ 1 + Gr/kj + w j 2 r 2 


(75) 


We note from Equation 75 ihat if the thermal time lag is zero or infinite, there is no effect on tlu* 
'‘steady state” natural frequency. For any finite value of r, the natural frequency is less than the 
frequency without the thermal effect. It can easily be shown that £2 has a minimum when 


WjT * 1 


(76) 


Using this expression in the more accurate Equation 72 results in 


2 W G 
n mia %w f 


2aj i In 


( 77 ) 


For the previously mentioned spacecraft parameters (except c and 5 } ) one can calculate the r and 6 } 
required for the observed morbit spin rate ripple period. A simple calculation using Equations 72 
and 76 to solve for G and hence $ s gives 

r * 2.04 seconds 
6 $ * 1.71 ft 

It is felt that the necessary time constant is rather short but may be realizable, while the steady 
state thermal extension is simply too large. The disparity between observed and calculated period 
could be due to several factors: this analysis has a very simple thermal model to evaluate thermal 
displacements and shadowing effects (both in a linear sense). A more elaborate thermal model 
may improve the comparison. Also, nonlinear dynamic coupling of the symmetric and antisym- 
metric modes could be occurring; or possibly uncertainties in mass properties and linear dimensions 
could be another. To explore the nonlinear effects due to dynamics and slightly different tip 
masses, Equations 21b, 24b, 30 and the thermal equations were solved on a digital computer. 

It was thought that perhaps the unobservable symmetric mode was large due to thermally driver, 
motion and through nonlinear coupling affect the period of oscillation of the smaller amplitude 
antisymmetric mode. Computer runs showed that for d { , 9, the order of 0. 16, 0.25 radians, 
respectively, the period of the antisymmetric mode was extended but the oscillations were no 
longer sinusoidal in form. Htis approach was dropped. 

A practical method to suppress thermally induced spin rate ripple is to utilize an appendage motion 
damper for the antisymmetric mode (Reference 1). It is a fluid filled ring, tuned to the antisym- 
metric mode frequency, similar to a passive r.utation damper but orientated so that the spin rate 
fluctuations excite fluid motion. This action dissipates energy through fluid viscous friction and 
eventually limits the spin rate ripple below discemability. This type of damper was flown on IMP-J, 
ISEE-A, C and DE-A. *11 spacecraft had long radial wire appendages subject to thermally induced 


» 1 * , 


spin rate rippie. There was no indication of spin rate ripple except going in and coming out of the 
earth’s shadow (ISEE-A and DE-A) which was rapidly eliminated by the damper. 

In summary, it is felt that the spin rate rippie anomaly observed on ISEE-B can be explained in 
terms of a thermally induced hypothesis that couples the ribbon angular motion with the thermally 
induced displacement. Analytic results show a possible exponential buildup of the inplane motion 
of the ribbon, ultimately limited by damping in the ribbon and spacecraft. The analysis indicates 
that the thermal lag time constant is fundamental for an increase in the antisymmetric mode period. 
Though there is good overall qualitative agreement with flight data, the numerical comparison of 
the ripple period is not as good as desired. 
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